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FIELD OF THE INVENTION 

The present invention relates to digital-analog converters and in particular to a 
method of generating a switching sequence for an unary array of conducting 
branches and a relative thermo metrically decoded digital-to-analog converter 
5 having a desired integral non-linearity (briefly INL) error function. 

Background of the invention 

Digital-analog (D/A) converters are widely used in many applications, such as 
wireless communications, signal reconstruction, waveforms generation etc.. 
Basically, a D/A converter may be realized with an array of conducting branches 
10 that generate binary weighted constant currents (binary array) and/or identical 
currents (unary array). These conducting branches may include capacitors or 
current generators, selectable by respective switches for delivering current toward 
a summing line or toward ground. 

For sake of simplicity, let us refer to a D/A converter, the conducting branches of 
15 which are realized with current generators, though what will be said will hold, 
mutatis mutandis, also for conducting branches realized with switched capacitors. 

The conducting branches of a binary array are selected by respective bits of an 
input bit string, while the conducting branches of an unary array are selected 
according to a fixed switching sequence for making the unary array generate a 
20 total current corresponding to the input digital value. 

For better comprehending the field of utility of this invention, let us refer to the 
so-called segmented converters, which generally have a binary weighted array of 
conducting branches, selected by the least significant bits of an input digital 
string, and an unary array that delivers a current corresponding to the digital value 
25 represented by the most significant bits, that is according to a so-called 
thermometer or thermometric decoding. 

A sample unary array is depicted in Figure 1. As may be noticed, the current 
generators are selected according to a certain switching sequence and the output 

l 



current lout corresponds to the sum of the currents circulating in the selected 
branches. Optionally, there may be a "dummy" area in the unary array, which may 
be used for realizing biasing circuits. In segmented digital-to-analog converters, 
the "dummy" area may be used for realizing a binary array. 

When the digital input value to be converted by the unary array represents a 
certain number n, the current generators of the unary array from 1 to n are 
switched on, thus generating an output current lout proportional to the digital 
input number n. 

Ideally, a D/A converter should generate an output signal that varies linearly with 
the input bit string, which would happen if all the conducting branches were 
identical. Unfortunately, mismatches between conducting branches due to 
inaccuracies of the fabrication process (process spread), make the current 
delivered by each branch not exactly equal to the design value, but affected by an 
error that may depend on the position of the conducting branch on the silicon 
substrate. 

In general, unary arrays are affected by a differential non-linearity (DNL) error 
and by an integral non-linearity (1NL) error. 

Indicating with 7 the average current delivered by the branches of an unary array 
and with Ij the current delivered by the>th branch of the array in the switching 
sequence, 

/ /= 7-(i + 0 

wherein e , is the relative deviation of the current Ij from the average current 7 . 

In an unary array not having a dummy area, the DNL error of the *-th branch in 
the switching sequence is 

; DNL(k) = e k 



2 



which represents a non-uniform deviation in the ideal current step amplitude 
between adjacent bit strings. 

The INL error function is defined as 

INLik^j^e; 

7=1 

and gives the deviation of the real analog output signal from its ideal value as a 
percentage of the average current 7 for any value of the switching sequence. 

The absolute INL error of a switching sequence is the maximum absolute value of 
the relative INL error function. 

An introduction on the INL and DNL errors of an unary array of a D/A 
thermometrically decoded converter is carried out in the article by Y. Cong and R. 
L. Geiger "Switching Sequence Optimization for Gradient Error Compensation in 
Thermometer-Decoded DAC Arrays", IEEE Trans, on circuits and systems - II: 
analog and digital signal processing, Vol. 47, No. 7, pages 585-595, July 2000. 

The DNL error of a branch can be reduced only by reducing the process spread 
and is independent from the switching sequence. By contrast, the INL error 
function strongly depends on the switching sequence, as it may be easily inferred 
from the following example. 

Two possible switching sequences of a mono-dimensional unary array having 
eight conducting branches are shown in Table 1. The conducting branch are 
affected by the indicated DNL errors (e). The absolute values of the underlined 
numbers are the absolute INL errors of the switching sequences. 

As may be noticed, the INL error function of the sequential switching sequence 
shows a maximum deviation of the analog output from its ideal value in the 
middle of the sequence. This situation is inconvenient, because the digital values 
input to a D/A converter are more likely in the middle of the range of conversion, 
rather than at the two ends thereof. Therefore, it is more convenient to switch the 
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mono-dimensional array of conducting branches of Table 1 according to a 
symmetrical sequence than according to a sequential sequence. Moreover, the 
absolute INL error for the symmetrical sequence is 7, while for the sequential 
sequence is 16. 

5 In this very simple case, it is possible to determine by successive trials the 
switching sequence with the smallest absolute INL error, but for two-dimensional 
unary arrays of about a thousand of conducting branches the number of 
combinations is too large for determining a switching sequence with the desired 
INL error function by trials. 

10 Many different methods of determining a switching sequence of a two- 
dimensional unary array of conducting branches and a relative D/A converter have 
been proposed. 

The patent US 6,118,398 by G.J. Fisher et al. discloses a digital-analog converter 
having an unary array of current sources that are selected according to a sequence 
15 that ensures a relatively small absolute integral non-linearity error. The suggested 
switching sequence is substantially a mixed symmetrical sequence, in which the 
current sources that are in the middle of the array have median positions in the 
switching sequence, while current sources that are in borderline regions of the 
array are at the beginning or at the end of the sequence. 

20 The patent US 5,057,838 by K. Tsuji et al. discloses a D/A converter having a 
plurality of conducting branches of a two-dimensional array, wherein the 
switching sequence is determined in order to make the center of the current 
contributions delivered by the conducting branches of the array coincide with the 
center of the array. 

25 The document "A 14-bit Intrinsic Accuracy Q2 Random Walk CMOS DAC" by 
Van der Plas, Steyaert et al., JSCC 12 Dec '99, discloses a method of determining 
the switching sequence of the switches of a D/A converter organized in matrix 
form, exploiting the so-called "Q 2 random walk" algorithm. 
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The document "Switching Sequence Optimization for Gradient Error 
Compensation in Thermometer-Decoded DAC Arrays", by Cong, Geiger, JSSC 7 
Jul '00, discloses an algorithm to find the so-called "INL bounded" switchmg 
sequence for an unary array of branches organized in matrix form. 
Unfortunately, the absolute INL errors of the D/A converters realized with the 
above techniques remain relatively large. 
OBJECT AND SUMMARY OF THE INVENTION 

As it will be more clearly described hereinafter, the array of conducting branches 
is affected by an error whose spatial distribution has an anti-symmetrical linear 
component, a quadratic symmetrical component, components of higher orders and 
a random component. The known thermometrically decoded converters are 
affected by quite large INL errors because the switching sequence of their unary 
array of conducting branches is not determined in a way that would optimally 
compensate the quadratic component of the error distribution. 

The objective of the method of this invention is to determine switching sequences 
of two-dimensional unary arrays of conducting branches of thermometrically 
decoded D/A converters, in a way that will ensure that the relative INL error 
function be contained between pre-established symmetrical upper and lower 
bound functions. These functions may be constant, such to ensure an absolute INL 
, error smaller than a pre-established value, or approaching zero in correspondence 
of mid way values of the switching sequence from a certain maximum value in 
correspondence of the two ends of the range of conversion. This last solution is 
preferable when the D/A converter is likely to work almost constantly about the 
middle of its range of conversion. 
5 When these upper and lower bound functions are constant, the obtained switching 
sequence compensates both the linear and the quadratic component of the error 
distribution and therefore is affected by a very small absolute INL error, which 
depends essentially on the random component of the error distribution. 
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The method of this invention may be easily implemented by a computer program 
and allows the realization of thermometrically decoded D/A converters affected 
by a known limited INL error function. 

The invention is more precisely defined in the annexed claims. 
Brief Description of the Drawings 

The advantages of this invention will become more evident through a detailed 
description referring to the attached drawings, wherein: 

Figure 1 illustrates a sample 16x16 unary array; 

Table 1 compares the INL error functions of two switching sequences for a 
mono-dimensional array; 

Figure 2 shows a sample joint error distribution; 

Figure 3 shows how to build a switching sequence according to the method of 
this invention; 

Figure 4 shows a flow chart of a preferred embodiment of the method of this 
invention; 

Figure 5 depicts a switching sequence determined by using the method of this 
invention for a 16x16 unary array; 

Figure 6 depicts a sequential switching sequence for a 16x16 unary array; 
Figure 7 depicts a symmetrical switching sequence for a 16x16 unary array; 
Figure 8 depicts a switching sequence for a 16x16 unary array of the 12-bit 
thermometrically decoded D/A converter embedded in the commercial devices 
MTC-xxl54, MTC-xxl74 and MTC-xx454 of Alcatel Microelectronics; 
Figure 9 depicts a switching sequence for a 16x16 unary array determined with a 
Q 2 random walk algorithm; 

Figure 10 depicts an "INL bounded" switching sequence for a 16x16 unary array; 
Figure 11 depicts a random switching sequence for a 16x16 unary array; 
Figure 12 depicts an anti-symmetrical switching sequence in respect to the center 
of a 16x16 unary array; 

Figure 13 depicts an improved anti-symmetrical switching sequence in respect to 



the center of a 1 6x 1 6 unary array; 

Figures 14 and 15 compare the WL error functions of the switching sequences of 

Figures from 5 to 13 in function of the direction 0of the linear error distribution; 

Tables 2 and 3 compare the maximum values of the absolute INL errors of many 
5 switching sequences obtained with known methods and the method of the 

invention for different values of the linear # and quadratic g q coefficients; 

Figure 16 depicts the switching sequence of an unary array of 1024 conducting 

branches of a 14-bit D/A converter of this invention, ten of which are 

thermometrically decoded. 
10 DESCRIPTION OF AN EMBODIMENT OF THE INVENTION 

Before illustrating the method of this invention, it is necessary to make some 
mathematical considerations. 

As stated hereinbefore, the conducting branches of an array are affected by an 
error, which in electronic equipments may be due to thermal and/or mechanical 
15 phenomena, doping differences, and other process spread mechanisms of 
inaccuracy, distributed according to an error distribution function over the array. 

In general, for a two-dimensional disposition of conducting branches, the relative 
error distribution e{x,y) may be approximated by a Taylor series expansion: 

s(x,y)=a 0 + a n x + a, 2 y + a 21 x 2 + a^y 2 + a 23 xy + --- 

20 Truncating this series at the second order, neglecting the term a^xy and 
supposing that a 2l =an, the linear ^(x,.y), quadratic e q {x,y) and joint £j(x,y) 
error distribution functions are, respectively, 

£ l (x,y)=g l (cos0 x + sm0 y); e q (x, y) = g q - (x 2 +y 2 )-a 0 ; 

e J (x,y)=e l (x,y)+e q (x,y) 
25 wherein 0 is the angle of the linear error gradient, g, is the slope, g q is a quadratic 
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coefficient and a 0 is an offset value. This joint error distribution, represented in 
Figure 2, has proven to be a good approximation of real error distribution 
functions of unary arrays. The difference between the real error distribution and 
the joint error distribution is due to the components of higher order thereof and to 
5 a random error component. 

In order to compensate errors due to the linear component of the error distribution, 
according to the method of the invention each pair of consecutive odd (2n-l) and 
even (2n) component of the switching sequence must be symmetrical in respect to 
the "center of gravity" of the two-dimensional array. In fact 

10 e l {x,y) = -e l (-x-y) 

and thus, if pairs of consecutive branches {a,b} and {g,h} are disposed as 
depicted in Figure 3, the linear component of the error introduced by each pair is 
null. 

Being 7~ q the continuous component of the quadratic error distribution over the 
15 area of the array of branches, that is the mean value thereof, the quadratic error 
distribution is also given by 

wherein e* c (x,y) is the alternate component of the quadratic error. 

The offset term T q does not contribute to the INL error because it causes only an 
20 offset error in the slope of the transfer function of the D/A converter. On the 
contrary, the alternate component has positive and negative values that make the 
transfer function nonlinear, thus originating the INL error. 

Therefore it is clear that, in order to have a switching sequence with a desired 
(small) INL error function, it is necessary to compensate the alternate component 
25 of the quadratic error distribution. 
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According to an essential embodiment of the method of the invention, the first 
step consists in defining an upper bound function and a lower bound function, 
symmetrical to each other, of the INL error function of the switching sequence to 
be generated. Then the error distribution function over the array is evaluated, in 
5 order to calculate the error associated to each pair of symmetrical branches. As 
stated before, each pair of successive branches must be symmetrical in respect to 
the center of gravity in order to compensate the component of the INL due to the 
linear error distribution. 

These bound functions define the range of variation of the INL error and they may 
10 be constant, if the absolute INL error must be minimized. 

The appropriate switching sequence is built by choosing a first pair of branches 
(1, 2) starting from any branch of the array and performing iteratively the 
following steps: 

- calculating a corresponding value of the INL error function of the switching 
15 sequence being built, 

- choosing as the successive pair, the pair of conducting branches that 
maximizes or minimizes the next value of the INL error function of the 
switching sequence though remaining comprised between the corresponding 
values of the upper bound and lower bound functions, 

20 - if all other pairs do not meet the preceding conditions, then repeating the steps 
restarting from the first by choosing every time a different first pair of 
branches, and if the conditions cannot yet be met, changing at least one of the 
bound functions and restarting from the first step. 

Finally, when the determination of the appropriate switching sequence has been 
25 completed, it is output. 

Preferably, the above described method is repeated for all possible first pairs of 
conducting branches, thus generating a set of switching sequences affected by an 
INL error function comprised between the same bound functions. The optimal 
switching sequence is chosen from this set according to a pre-established 
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criterion. 



Optionally, the bound functions may be closer to zero in correspondence of mid 
way values of the range of the converter compared to their value in 
correspondence of the two ends of the range of the converter. 

5 A preferred embodiment of the method of this invention is described by the flow 
chart of Figure 4. In this case the symmetrical bound functions are constant and 
this constant value is half the maximum absolute value of the DNL error of the 
array. 

An example of a 16x16 matrix describing a switching sequence determined 
10 according to the method of the invention illustrated in Figure 4 and for gi = g q = 
0.5 is depicted in Figure 5. As may be noticed, the switching sequence has been 
determined by selecting every even conducting branch (2, 4, 6 ...) symmetrical to 
the preceding odd branch (1, 3, 5 .. .) in respect to the center of the array. 

Preferably, but not necessarily, the map of the switching sequence will be anti- 
15 symmetrical in respect to an axis of symmetry of the array, that means that the last 
branch of the sequence is symmetrical to the first branch in respect to that axis of 
symmetry, and so forth for the other branches. For example, the positions of the 
branches 256, 255, 254 etc. of the switching sequence of Figure 5 are symmetrical 
to the positions of the branches 1, 2, 3 etc. in respect to the horizontal axis of 
20 symmetry of the array. 

For comparison purposes, eight other matrices obtained with different algorithms 
are depicted in Figures from 6 to 13. 

The resulting absolute INL error values in function of the angle 0 of the gradient 
of the linear error are compared in Figures 14 and 15 for gi = g q = 0.5. The values 
25 have been calculated after having normalized to 1 the maximum value of the 
linear or quadratic errors. As may be noticed, the absolute INL error of the 
switching sequence of Figure 5 is always smaller than that of the other sequences. 
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Moreover, differently from the switching sequences of Figures from 6 to 13 
obtained by using the known methods, the absolute INL error of the switching 
sequence of Figure 5 is substantially independent from the angle 9. This is very 
important because it ensures that the method of the invention is not tied to a 
particular shape of the error distribution, as the method disclosed in the 
aforementioned paper by Y. Cong and R. L. Geiger. 

Table 2 compares many switching sequence, indicating the maximum value of the 
absolute INL error for each of them. The references [1] and [2] indicate that the 
switching sequence is obtained using the method described in the mentioned 
paper by G. A. Van der Plas et al. and by Y. Cong and R. L. Geiger, respectively. 

As may be noticed, the switching sequence of Figure 5 is the best one. This result 
is confirmed even using different values of gi and g q , as shown in Table 3. Even if 
these parameters undergo sensible variations, the maximum absolute INL error of 
the switching sequence obtained with the present method is substantially 
independent from them. This extraordinary result confirms that the method of this, 
invention effectively compensates both the linear and the quadratic component of 
the error distribution. *■ 

The method of this invention may be applied whichever the shape of the array of 
conducting branches is, not only to two-dimensional square arrays. 

A switching sequence generated according to the preferred embodiment of the 
method illustrated in Figure 4 for a substantially oval unary array of conducting 
branches of a 14-bit D/A converter of the invention, ten of which are 
thermometrically decoded, is depicted in Figure 16. As may be noticed, the 
switching sequence is anti-symmetrical in respect to the vertical axis of symmetry. 
The ten most significant bits of the input bit string of the D/A converter of the 
invention are thermometrically decoded, while the four least significant bits select 
as many conducting branches of a binary scaled array. 
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CLAIMS 

1. A method of generating a switching sequence for an unary array of 
conducting branches of a thermo metrically decoded digital-to-analog converter, 
the conducting branches of which are affected by an error distributed over the 
array according to a certain error distribution function, comprising the steps of: 

preliminary evaluating said error distribution function, fixing an upper 
bound function and a lower bound function symmetrical to each other 
of the INL error function associated to the switching sequence, 
determining the center of gravity of the array and calculating, by said 
preliminarily evaluated error distribution function, error values 
associated to each pair of conducting branches symmetrical in respect 
to the center of gravity of the array; 
choosing pairs of successive conducting branches in the switching sequence 
such that every even conducting branch (2n) is substantially 
symmetrical to the preceding odd conducting branch (2n-l) in respect 
to said center of gravity; 
building a switching sequence with an INL error function comprised 

between said upper bound and lower bound functions by: 
a) starting from a first chosen pair of conducting branches (1, 2) of the 
switching sequence iteratively performing the following steps from bl) 
to b3) for all other pairs of conducting branches of the array: 
bl) calculating a corresponding value of the INL error function of the 

switching sequence being built, 
b2) choosing as the successive pair, the pair of conducting branches that 
maximizes or minimizes the next value of the INL error function of the 
switching sequence though remaining comprised between the 
corresponding values of said upper bound and lower bound functions, 
b3) if all said other pairs do not meet the conditions of point b2), then 
repeating the steps from a) to b3) choosing every time a different first 
pair of branches, and if the conditions cannot yet be met, changing at 



least one of said bound functions and restarting from point a); 
c) when all other pairs of conducting branches of the unary array meet the 
conditions, outputting the resulting switching sequence. 

2. The method of claim 1 , further comprising the operations of: 
5 repeating the steps from a) to c) for all possible first pairs of conducting 

branches using the same bound functions, obtaining a set of switching 
sequences whose INL error function is comprised between said bound 
functions; 

choosing the switching sequence of said set that satisfies a pre-established 
10 criterion. 

3 The method of claim 1 , wherein said bound functions are constant. 

4. The method of claim 3, wherein the value of said upper bound function 
is half the maximum value of the absolute value of the DNL error affecting the 
conducting branches of the array. 
15 5. The method of claim 1, wherein said bound functions are closer to 0 in 

correspondence of mid way values of said switching sequence than in 
correspondence of the two ends of a range of conversion. 

6. The method of claim 2, wherein said criterion consists in choosing the 
switching sequence of said set that is affected by the smallest absolute INL error. 

20 7. The method of claim 2, wherein said criterion consists in choosing the 

switching sequence of said set that is affected by the smallest absolute DNL error 
in correspondence of raid way values of the switching sequence. 

8. The method of claim 2, wherein said criterion consists in choosing the 
switching sequence of said set that is affected by the most evenly oscillating INL 

25 error function. 

9. The method of claim 1, wherein said pairs of conducting branches are 
also chosen to make the switching sequence being built anti-symmetrical in 
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respect to an axis of symmetry of the unary array. 

10. A thermometrically decoded digital/analog converter, comprising an 
unary array of conducting branches selectable by respective switches, 
characterized in that the switching sequence of said branches is generated by 
using the method of claim 1 . 

11. A computer program loadable in an internal memory of a computer, 
comprising a software code for performing the steps of the method of claim 1 
when said program is executed on a computer. 
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"METHOD OF GENERATING A SWITCHING SEQUENCE FOR AN 
UNARY ARRAY OF CONDUCTING BRANCHES AND A RELATIVE 
THERMOMETRICALLY DECODED DIGITAL-TO-ANALOG 

CONVERTER" 

ABSTRACT 

A method for determining switching sequences of two-dimensional unary arrays 
of conducting branches of thermometrically decoded D/A converters, in a way 
that will ensure that the relative INL error function be contained between pre- 
established symmetrical upper and lower bound functions, has been found. When 
these, upper and lower bound functions are constant, the obtained switching 
sequence compensates both the linear and the quadratic component of the error 
distribution and therefore is affected by a very small absolute INL error, which 
depends essentially on the random component of the error distribution. 

This method may be easily implemented by a computer program and allows the 
realization of thermometrically decoded D/A converters affected by a known , 
limited INL error function. 
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Fig. 15 
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16 x 16 matrix comparison result (gp g Q =0.5) 



Row-Col 



Hierarchic. 



Improvements 



Algorithm/Sequence 
Row-Column basic 
Row-Column symm. Seq. 
Row-Column hier symm seq. [2] 
12-bit C05 DAC 

Q2 random walk [1] 

Hier. switch, IN L-bounded [2] 
Random (1 trial) 
Random (best of 1k trials) 
Magic 

Modified Magic 

This work 

TAB. 2 



Maximum absolute INL 
38.1 
20 
10.2 
5.9 

2.8 

1.93 
6.8 
3.4 
10.2 
2 6 
1).96 



^ Maximum absolute INL error 


9i = 


0.3 


0.5 


0.7 


9q = 


0.7 


0.5 


0.3 


Q z rand walk 


2.26 


2.8 


3.38 


INL-bounded 


1.56 


1.93 


2.4 


This work 


0.95 


0.96 


1.04 



Tab- 3 
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